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NOTATION 

A  The  surface  area  of  a  frustum,  ft^ 

C  Maximum  radius  of  a  body,  ft,  C  »  d/2 

CD  Drag  coefficient 

CH  Head  coefficient,  H/q0 

Cq  Mass  flow  coefficient 

Cp  Pressure  coefficient 

ft 

Cp  Ideal  forebody  power  coefficient 

Cv  Prismatic  coefficient,  fraction  of  usable  volume  contain  by 

47 

the  body,  Cv  -  - - - 

ird2*. 

d  Maximum  diameter  of  body,  ft 

H  Fluid  total  head,  lb/ft2  gage 

Kj  Maximum  velocity  ratio 

K2  Final  velocity  ratio 

L , i  Body  length,  ft 

J/d  Fineness  ratio 

M  Fraction  of  body  length  where  body  is  thickest 

n  Exponent  used  in  velocity  curve  fit 

P  Power,  ft-lb/sec 

Ps  Static  pressure 

q  Dynamic  pressure,  lb/ft4,  q  ■  1/2P  'J 

RQ  Radius  of  curvature  at  the  nose,  ft 

rQ  Nondimensional  radius  of  curvature  at  the  nose,  rQ  ■  RQ  •  i/d 

r|  Body  radius  at  the  location  of  the  boundary  layer  rake,  ft 

Re ^  Reynolds  number  based  on  volume  raised  to  the  one  third  power 
S  Surface  distance,  ft 
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V  Velocity,  ft/s er 

y  Height,  ft 

5  Boundary  layer  height,  ft 

V  Volume,  ft  ,  (for  BLC  bodies,  the  portion  of  the  body  volume  upstream 

of  the  inlet  location;  otherwise  total  body  volume) 

6  Angle  between  normal  to  body  surface  and  the  vertical,  degrees 

p  Density,  slugs/ft-* 

v  Kinematic  viscosity,  ft^/sec 

Subscripts 

i  Value  at  heights  normal  to  body,  y^^ 

1  Local  conditions  at  station  1,  the  assumed  inlet  location 

0  Free-stream  conditions 

te  Trailing  edge  of  the  body 
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ABSTRACT 


One  objective  of  the  Advanced  Undersea  Vehicle  (AUV)  program 
is  to  design  a  low  drag  vehicle*  The  approach  to  this  objective 
in  this  investigation  is  boundary  layer  control  by  means  of  an 
annular  suction  slot  located  on  the  afterbody.  Although  wind 
tunnel  data  showed  significant  reduction  in  propulsive  power 
over  conventional  shapes,  an  attempt  was  made  to  achieve  further 
reduction  by  means  of  forebody  shaping.  Two  methods  were  used 
to  vary  the  geometric  parameters  for  this  analysis.  The  direct 
method,  based  on  the  mathematical  development  of  the  Series  58 
bodies,  allows  the  definition  of  a  shape  by  a  fifth-order 
polynomial  based  on  the  four  fundamental  parameters  of  fineness 
ratio,  nose  radius  of  curvature,  location  of  maximum  thickness, 
and  prismatic  coefficient.  The  inverse  method  allows  various 
velocity  distributors  to  define  the  body  shape.  The  shapes 
derived  by  this  method  have  flat  velocity  distributions  and 
show  similar  trends  as  the  polynomial  shapes  (about  3-percent 
reduction  in  propulsive  power).  The  range  of  fineness  ratios 
analyzed  was  from  i  to  10  at  a  volume-based  Reynolds  number  of 
3.2  million.  In  the  range  of  2.5  to  8,  fineness  ratio  did  not 
affect  propulsive  power  more  than  6  percent.  A  maximum  improve¬ 
ment  of  3  percent  was  shown  by  varying  the  meridian  section. 


ADMINISTRATIVE  INFORMATION 

This  work  was  performed  by  the  New  Vehicle  Office  (Code  1603)  of  the 
Aviation  and  Surface  Effects  Department  at  the  David  W.  Taylor  Naval  Ship  Research 
and  Development  Center  (DTNSRDC).  The  project  was  funded  by  the  Defense  Advanced 
Research  Projects  Agency  (DARPA)  under  Program  Element  62702E,  DTNSRDC  Work  Unit 
1603-604. 


INTRODUCTION 

One  aspect  of  the  Advanced  Undersea  Vehicle  (AUV)  program  addresses  the 
development  of  low-power  requirements  for  an  axisymmetric  turbulent  flow  hydro- 
dynamic  body.  The  present  concept  uses  the  ingestion  of  the  boundary  layer  fluid 
and  the  reenergization  of  the  fluid  to  propel  the  vehicle.  By  ingesting  the  fluid 
at  the  proper  location  and  with  proper  afterbody  shaping,  flow  separation  can  be 
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avoided  even  with  very  full  afterbodies.  The  ability  to  design  full  afterbodies 
without  a  power  penalty  has  resulted  in  high  prismatic  coefficient  bodies 
(Cy  >  0.8)  with  low-power  requirements.  Analytical  studies  have  been  performed  to 
develop  body  shapes  for  experimental  verification  in  a  wind  tunnel.**  The  body 
shape  selected  for  thl9  experiment  uses  a  Reichardt  nose  and  parallel  mid-bodies 
which  varied  the  body  length-to-ciameter  ratio  (4,  5,  and  6);  see  Figure  1.  The 
full  afterbody  shape  was  selected  based  on  a  preliminary  design  for  the  prupulsor 
and  shroud . 

The  three  designs  were  experimentally  evaluated  in  a  wind  tunnel  in  an 
unpowered  condition  using  a  fan  mounted  downstream  of  the  model  connected  by  a  duct 
to  provide  the  suction  air  simulating  the  propulsor.  Pressure  measurements 
provided  data  to  determine  the  power  coefficients  for  the  bodies  at  various  suction 
flow  rates.  One  of  the  many  observations  during  the  wind  tunnel  investigation  was 
that  the  body  with  a  fineness  ratio  6  (l.e.,  length-to-maximum-diameter  ratio) 
performed  better  than  the  lower  fineness  ratio  bodies.  This  study  was  undertaken 
to  verify  the  wind  tunnel  results  and  to  determine  if  another  shape  could  perform 
even  better  than  the  shapes  selected  for  the  wind  tunnel  experiments. 

One  way  of  assessing  the  performance  of  a  vehicle  is  by  calculating  the  power 
required  to  propel  the  vehicle  at  a  constant  speed,  that  is,  propulsive  power.  An 
attempt  was  made  to  reduce  propulsive  power  on  the  wind  tunnel  model  by  means  of 
forebody  shaping.  It  was  believed  that  an  efficient  means  of  accelerating  the 
fluid  over  the  forebody  such  that  boundary  layer  growth  would  be  controlled  would 
result  in  a  reduction  in  required  power.  Forebody  shaping  to  reduce  drag  had  been 
investigated  a  number  of  times,  but  the  feasibility  of  efficient  afterbody  boundary 

*A  complete  Hating  of  references  is  given  on  page  33. 
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EL  CEILING 


Figure  1  -  Component  Schematic  of  Boundary  Layer  Control  Model 

Installed  in  Wind  Tunnel 


't’RMnp 


layer  control,  which  was  demonstrated  earlier  in  the  program,  provides  two  new 
considerations  that  suggest  a  separate  study.  First,  because  separation  is 
avoided,  i.e.,  no  form  drag,  the  primary  component  of  drag  to  be  considered  is  skin 
friction.  In  addition,  since  an  afterbody  is  required  for  boundary  layer  control, 
the  analysis  of  the  skin  friction  should  be  conducted  only  from  the  nose  of  the 
body  to  a  ooint  where  an  inlet  would  be  placed.  Second,  only  turbulent  flow  is  of 
interest  in  this  study  whereas  some  other  studies  dealt  with  laminar  flow.  These 
tvTO  issues  justified  proceeding  with  a  forebody  shaping  program  at  DTNSRDC.  (See 
References  3  and  4  for  details  of  other  studies.) 

EVALUATION  CRITERIA 

It  is  assumed  that  with  a  properly  developed  BLC  aftbody  the  form  drag  will 
be  negligible,  external  aftbody  friction  will  be  very  small  and  its  dependence 
on  forebody  geometry  can  be  neglected,  and  the  internal  efficiency  with  which 
boundary  layer  fluid  can  be  ingested,  reenergized  and  exhausted  is  also  indepen¬ 
dent  of  forebody  geometry.  These  assumptions  greatly  simplify  the  problem  of 
comparing  the  merits  of  different  overall  body  geometries  since  they  can  be  com¬ 
pared  on  the  basis  of  boundary  layer  properties  at.  an  assumed  propulsor  Inlet 
location  and  the  details  of  aftbody  internal  and  external  shaping  need  not  be 
considered. 

The  inlet  location  was  selected  as  a  point  slightly  upstream  from  the  occur¬ 
rence  of  boundary  layer  separation  as  determined  from  a  calculation  of  the  bound¬ 
ary  layer  development  under  the  influence  of  the  calculated  potential-flow  pressure 
distribution.  By  choosing  to  nondimenslonalize  results  in  terms  of  the  body  vol¬ 
ume  upstream  of  the  chosen  inlet  location,  it  was  found  that  the  effect  of  small 
variations  in  the  exact  choice  of  inlet  location  was  made  negligible.  (This  choice 
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of  reference  volume  was  necessary  In  any  case  in  order  to  avoid  having  to  deal 
with  the  details  of  aftbody  geometry.) 

The  quantity  chosen  as  the  measure  of  relative  merit  of  the  various  body 
geometries  is  ideal  forebody  power  coefficient ,  Cp*.  a  nondimensional  measure  of 
the  power  required  by  an  ideal  propulsor  to  restore  the  boundary  layer  fluid  at 
the  chosen  inlet  location  to  free-stream  total  head. 

This  coefficient  is  defined  as: 


where 


* 


P 

«oV 


2/3 


6 

/  Vi  (%-Hi>  dAi 
o 


(1) 

(2) 


dAj^  *  2tt  (rj  +  y^  cos  0)  dy 


(3) 


But 


Vi  dAi  -  V2/3  dCQ 


(4) 


Hence 


)  dCq 


(5) 


The  requisite  boundary  layer  information,  namely  CH  as  a  function  of  Cq. 
was  obtained  from  the  Transition  Analysis  Program  System  (TAPS)^’^  which  was  also 
one  of  the  primary  tools  used  in  the  analysis  and  development  of  the  earlier 
wind  tunnel  configurations.  Correlation  between  TAPS  and  actual  wind  tunnel 
data  has ,  to  some  degree,  been  completed.  Very  good  comparisons  between  the 
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TAPS  potential  solution  and  pressure  coefficient,  (p  -  p  )/q  ,  measurements 

"1  8o  ° 

have  been  found.  In  Figure  2,  a  sample  plot  shows  the  pressure  distribution 
over  the  wind  tunnel  model  (  i/d  -  6)  with  the  TAPS  prediction.  Also,  good  agree¬ 
ment  has  been  found  between  the  predicted  energy  content  of  the  fluid  just  up¬ 
stream  of  the  inlet  and  the  measured  energy  (head  profile).  Predicted  head 
profile  is  compared  with  the  measured  head  profile  in  Figure  3.  The  head  pro¬ 
file  is  presented  by  a  head  coefficient  versus  a  mass  flow  coefficient;  the 
latter  is  similar  to  the  stream  function  concept,  that  is,  the  mass  flow  is 
constant  between  streamlines.  This  method  of  presentation  is  essential  because 
the  experimental  data  revealed  significant  variation  of  static  pressure  across 
the  boundary  layer.  Since  TAPS  assumes  uniform  static  pressure  across  the  bound¬ 
ary  layer,  comparisons  of  H  vs  y  and  V  vs  y  separately  would  suggest  very  poor 
correlation  even  though  the  actual  correlation  in  terms  of  flow  quantities  needed 
for  this  study  is  excellent. 

The  wind  tunnel  model  geometries  were  also  analyzed  by  TAPS  at  the  same  con¬ 
ditions  as  the  wind  tunnel  data  points  so  that  further  correlation  could  be  made 
between  predictions  and  actual  measurements.  The  wind  tunnel  models  had  a  13-tube 
rake  of  which  11  tubes  were  for  total  pressure  and  2  were  for  static  pressure 
measurements  (Figure  4).  This  rake,  located  just  forward  of  the  inlet,  is  referred 
♦•o  as  the  boundary  layer  rake.  The  boundary  layer  rake  measurements  were  integra¬ 
ted,  as  shown  in  Equations  (1)  through  (5),  in  order  to  correlate  the  analysis  with 
actual  data.  The  results  of  this  correlation  showed  very  good  agreement,  as  indi¬ 
cated  in  Table  1 . 


GEOMETRY  DEVELOPMENT 

In  order  to  compare  the  different  shapes  and  to  understand  their  variation  in 
performance,  a  method  was  sought  which  would  allow  the  variation  of  geometric 
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Distribution  Over 


Figure  3  -  Energy  Content  versus  Mass  Flow  in  Boundary  Layer 


Figure  4  -  Boundary  Layer  Rake  Schematic 
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TABLE  1  -  CORRELATION  BETWEEN  TAPS  PREDICTION  AND  WIND  TUNNEL  DATA 


Wind  Tunnel 
Model,  t/d 

Predicted 

V 

Measured 

V 

Variation 
( X ) 

6 

0.0201 

0.0206 

2.5 

5 

0.0202 

0.0200 

1.1 

4 

0.0195 

0.0193 

1.0 

parameters.  These  parameters  include  fineness  ratio,  prismatic  coefficient,  nose 
radius  of  curvature,  and  the  location  of  maximum  thickness.  Two  approaches  are 
available  to  accomplish  this.  The  first  approach,  the  direct  method,  is  to  define 
a  geometry  based  on  the  four  parameters  and  then  analyze  its  flow  characteristic. 

The  second  approach,  the  inverse  method,  is  to  prescribe  a  velocity  distribution 
and  find  the  geometry  that  matches. 

DIRECT  METHOD 

The  direct  method  is  a  modification  to  the  work  done  by  Gertler  and  Landweber^ 
in  1950.  The  problem  simply  stated:  A  polynomial  is  desired  which  defines  the 
geometry  of  a  body  of  revolution.  It  is  desired  that  four  parameters  be  available 
for  the  user  to  vary,  thereby  generating  a  family  of  geometrically  similar  bodies. 
The  parameters  are  fineness  ratio,  prismatic  coefficient,  nose  radius  for  curvature, 
and  the  location  of  maximum  thickness.  The  equation  takes  the  form: 

y2(x)  -  2  rD  R(x)  +  C2Q(x)  +  CyP(x)  (6) 

where  R(x),  Q(x),  and  P(x)  are  fifth-order  polynomials. 

The  one  parameter  that  is  not  apparent,  but  of  major  importance,  is  the 
location  of  maximum  thickness  of  the  body,  M.  Due  to  the  mathematical  limitations 
of  Equation  (6),  M  was  chosen  to  be  0.7  of  the  body  length  for  all  geometries 
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analyzed.  The  derivation  of  R(x),  Q(x),  and  P(x)  la  shown  in  the  appendix  and  is 
similar  to  the  method  presented  in  Reference  7. 

The  major  portion  of  this  study  was  done  with  geometries  generated  by  the 
direct  method.  Certain  restrictions  are  imposed  in  Equation  (6)  in  order  to 
ensure  a  smooth  body.  An  example  of  these  limitations  is  shown  in  Figures  5  and 
6.  Figure  5  shows  a  region  within  which  only  one  point  of  zero  slope  exists  on  the 
body  contour.  This  alone  does  not  ensure  a  realistic  geometry  nor  a  smooth  one. 
Figure  6  shows  regions  in  which  y  <  0  or  y"  >0,  which  would  indicate  a  nonrealis- 
tic  geometry  and  a  geometry  with  inflection  points,  respectively.  An  example  of 
a  polynomial  generated  body  is  shown  in  Figure  7  for  a  fineness  ratio  5  body. 

For  a  fineness  ratio  of  4,  a  range  of  prismatic  coefficients  and  radii  of 
curvature  were  tested.  The  investigation  of  fineness  ratio  4  geometries  was  done 
to  understand  the  influence  that  the  prismatic  coefficient  and  nose  radius  of 
curvature  have  on  the  power  required  to  self-propel  the  bodies. 

The  very  nature  of  this  method  enables  it  to  be  converted  easily  to  a 
computer  program.  After  conversion,  a  systematic  procedure  was  developed  to  start 
with  a  geometry  and  result  in  a  final  Cp*.  This  method  was  also  extended  from  a 
fifth-order  to  a  sixth-order  polynomial.  As  no  significant  Improvement  was  shown, 
this  was  not  pursued. 

INVERSE  METHOD 

Included  in  the  TAPS  package  is  an  inverse  potential  flow  program  developed 
by  James®  and  the  McDonnell-Douglas  Corporation.^’*®  Input  to  the  program  is  a 
velocity  distribution,  and  by  means  of  iteration  a  resulting  geometry  is  output. 
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Figure  7  -  Polynomial  Generated  Geometry 


The  following  are  sample  input: 


where  S  ■  1  at  the  location  of  the  inlet,  and  S  ■  0  at  the  nose. 

By  varying  K^,  K2»  n»  and  Ste  a  family  of  velocity  distributions  can  be 
generated  (Figure  8). 


SURFACE  DISTANCE,  S 

Figure  8  -  Input  Velocity  Distribution 


A  number  of  attempts  were  made  using  the  TAPS  inverse  code,  but  the  approach 
was  discontinued  due  to  unsatisfactory  results.  Subsequent  analysis  of  an  output 
geometry  did  not  compare  with  the  input  velocity  distribution.  The  failure  of  the 
TAPS  inverse  method  is  believed  due  to  the  convergence  logic  and  its  particular 
coding;  therefore,  a  second  inverse  method  was  developed.  This  method,  based  on  a 
surface  discrete  vortex  distribution,  was  originally  coded  on  a  Tektronix  terminal 
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for  two-dimensional  flow.  Later,  the  method  was  developed  for  an  axisymmetric  flow 
and  finally  coded  in  FORTRAN  for  use  on  the  CDC-6600-type  computers . 

The  program  receives  as  input  a  function  for  the  body  length.  The  function 
used  was  a  series  of  straight  lines;  see  Figure  9.  The  broken  line  in  the  figure 
results  in  a  cusp-type  of  afterbody  and  is  similar  in  shape  to  the  velocity 
distribution  for  such  a  body.  The  solid  line  represents  a  symmetric  velocity 
distribution  (the  type  of  input  used  in  this  analysis).  The  resulting  geometry 
is  symmetrical,  which  would  change  in  shape  based  on  changes  in  the  slopes  in 
the  fore  and  aft  input  distribution.  The  resulting  geometries  were  checked  in  the 
TAPS  potential  flow  code  to  verify  a  flat  velocity  distribution.  A  sample  result 
is  shown  in  Figure  10. 


Figure  9  -  Typical  Input  Function  for  Inverse  Code 


A  number  of  other  geometries  were  analyzed  to  put  the  final  results  in 
perspective.  The  laminar  flow  forebody  shape  of  Goldschmied1 1  was  analyzed  along 
with  the  wind  tunnel  model  of  fineness  ratios  of  4,  5,  and  6. 
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RESULTS  AND  DISCUSSION 


Ac  low  fineness  ratios,  the  component  of  drag  due  to  flow  separation  (residual 
drag)  becomes  more  important.  This  phenomenon  is  seen  in  Figure  11  where  the 
total  drag  of  a  typical  Series  58  body  is  shown  along  with  the  components  of  drag, 
skin  friction,  and  residual.  Note  that  at  the  low  fineness  ratios  the  residual  drag 
controls  the  shape  of  the  total  drag  curve.  By  applying  afterbody  boundary  layer 
control,  the  residual  drag  component  can  be  eliminated,  and  the  optimum  fineness 
ratio  for  reducing  propulsive  power  can  be  found  based  on  a  tradeoff  of  skin 
friction  drag  alone. 

When  attempting  to  reduce  the  skin  friction  of  a  body  of  revolution,  there  is 
a  distinct  tradeoff.  By  reducing  fineness  ratio,  the  wetted  area  required  to 
enclose  a  given  body  is  reduced  as  shown  in  Figure  12.  This  should  reduce  skin 
friction  drag.  However,  reducing  fineness  ratio  also  increases  the  flow  velocities 
over  the  surface,  increasing  skin  friction  drag.  An  analysis  was  conducted,  using 
the  methods  previously  described,  to  define  the  fineness  ratio  which  optimizes  the 
reduction  in  wetted  area  with  the  increase  in  surface  velocity.  (An  assessment  of  a 
low  drag  spherical  submersible  was  conducted  at  DTNSRDC  in  1981.*) 

The  analysis  was  conducted  as  a  two-step  process.  First,  for  a  single 
fineness  ratio,  (i/d  ■  4),  the  influence  of  prismatic  coefficient  and  nose  radius 
of  curvature  was  determined.  The  results  in  Table  2  show  a  range  of  less  than  2 
percent.  Despite  such  a  small  range,  it  was  apparent  that  the  higher  prismatic 
coefficient  and  nose  radius  of  curvature  bodies  resulted  in  the  lower  power 
coefficients.  This  was  the  assumption  used  for  the  other  fineness  ratios. 


*As  documented  by  Leitner  and  McCabe  in  a  report  of  higher  classification. 
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_ ELLIPSOID  WETTED  AREA _ 

WETTED  AREA  OF  SPHERE  OF  EQUAL  VOLUME 


TABLE  2  -  PROPULSIVE  POWER  AS  A  FUNCTION  OF  PRISMATIC  COEFFICIENT 
AND  NOSE  RADIUS  OF  CURVATURE 
U/d  -  A) 


CP* 

CV 

ro 

0.01907 

0.55 

0.005 

0.01897 

0.60 

0.005 

0.01893 

0.015 

0.01894 

0.65 

0.005 

0.01887 

0.015 

0.01883 

0.025 

0.01897 

0.70 

0.005 

0.01893 

0.015 

0.01886 

0.020 

0.01879 

0.035 

0.01884 

0.75 

0.030 

0.01882 

0.035 

The  second  step  was  to  determine  the  influence  of  fineness  ratio.  The  results 
of  this  analysis  are  shown  in  Figure  13  and  Table  3.  The  primary  method  used  was 
the  polynomial  generated  shapes  which  make  up  the  base-line  curve.  The  polynomial 
shapes  range  from  fineness  ratios  of  1  up  to  10.  Figure  13  also  includes  the 
various  geometries  used  in  the  program  analysis.  The  wind  tunnel  model  with  the 
Reichardt  nose  displays  somewhat  lower  values  than  the  base  line.  A  third  family 
of  geometric  shapes  is  shown  with  still  lower  propulsive  power  coefficient.  These 
shapes,  representing  a  flat  pressure  distribution,  show  an  improvement  of  approxi¬ 
mately  3  percent  from  the  base  line  with  a  smaller  improvement  over  the  wind 
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POWER  COEFFICIENT,  C 
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tunnel  model.  The  characteristic  pressure  distributions  of  these  four  geometries 
are  shown  in  Figure  14. 

TABLE  3  -  TABULATED  RESULTS  OF  POWER  COEFFICIENT  VERSUS  FINENESS  RATIO 


i/d 

S* 

Geometry  Type 

1 

0.02765 

Polynomial  geometry 

2 

0.02137 

3 

0.01920 

4 

0.01880 

5 

0.01880 

b 

0.01910 

7 

0.01955 

8 

0.01990 

10 

0.02055 

4 

0.01850 

Wind  tunnel  model 

5 

0.01860 

b 

0.01895 

2.8 

0.01870 

Symmetric  body 

3.5 

0.01830 

4.9 

0.01835 

3 

0.01950 

Goldschmied 

CONCLUSIONS 


Although  the  optimum  fineness  ratio  appears  to  be  4,  the  curve  is  quite  flat. 
An  important  consequence  of  this  is  that  other  design  issues,  which  may  necessitate 
changes  in  the  vehicle's  length  or  diameter,  can  be  accommodated  with  impunity  over 
a  range  of  fineness  ratios  (2.5  to  8).  (In  practical  terms,  this  relative 
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Figure  1A  -  Pressure  Distribution  of  Various 
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insensitivlCy  might  extend  to  even  lower  fineness  ratio  bodies  since  the  imposition 
of  transition  to  turbulent  flow  at  10  percent  body  length  may  be  overconservative 
for  them.  The  sphere,  for  example,  has  a  powerful  favorable  pressure  gradient 
extending  much  further  than  10  percent  of  length  which  might  maintain  a  signi¬ 
ficantly  greater  region  of  laminar  flow  even  in  rather  unfavorable  practical 
conditions. ) 

For  a  constant  fineness  ratio,  the  different  shapes  varied  by  at  most  about 
4  percent  of  Cp*.  Relative  to  the  Reichardt  shape  selected  for  the  wind  tunnel 
model  improvements  of  1.5  percent  were  obtained. 

By  application  of  boundary  layer  control  to  underwater  submersibles ,  a  low- 
power  vehicle  can  be  designed  with  a  fineness  ratio  of  2.5  and  still  have  a  high 
prismatic  coefficient. 
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APPENDIX 

FIFTH-ORDER  POLYNOMIAL  DERIVATION 

The  form  of  equation  y2(x)  ■  2rQR(x)  +  CyP(x)  +  C2Q(x),  which  is  the 
sectional  area  curve,  is  of  the  fifth  order.  The  resultant  meridian  section 
can  be  easily  defined  hy  taking  the  square  root  of  y2.  Substituting,  f(x)  fur 
y2(x)  gives  the  following  boundary  conditions: 

a.  f(x  ■  0)  -  0 

b.  f(x  -  1)  -  0 
df 

c.  —  (x  -  0)  -  2rQ 
dx 

d.  f  (x  -  M)  -  C2 
df 

e.  —  (x  ■  M)  ■  0 
dx 

1 

f.  /  f(x)  dx  -  C2  Cv 
o 

where  M  -  location  of  maximum  thickness. 

Since  f(0)  -  0,  regardless  of  the  values  of  rQ,  C,  and  Cy, 

R( 0)  -  P(0)  -  Q( 0)  «  0  (1) 

Similarly,  Equations  (1)  through  (6)  are  obtained  from: 

Condition  b 

R( 1 )  -  P(l)  -  0(1)  -  0  (2) 

Condition  c 

R'(0)  -  1,  P  '(0)  -  Q'(0)  -  0  (3) 

where  the  prime  denotes  differentiation  with  respect  to  x. 

Condition  d 

R(M)  -  P(M)  -  0,  Q(M)  -  1  (A) 

Condition  e 

R'(M)  -  P'(M)  -  Q'(M)  -  0  (5) 
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Condition  f 

1  1  1 

f R(x)dx  -  f  Q(x)dx  -  0,  J  P(x)dx  -  C2 
0  0  0 


(6) 


The  values  of  R(x),  P(x),  and  Q(x)  are  derived  on  the  basis  of  the  relation  In 
Equations  (1)  through  (6). 


EVALUATION  OF  R(x) 

Since  R(0)  •  R(M)  *  R*"(M)  ■  R(l)  ■  0,  and  R(x)  is  a  polynomial  of  the  fifth 
degree,  Equation  (7)  becomes 

R(x)  -  x(x  -  1)  (x-M)2(a0  +  ottx)  (7) 

The  terms  aQ  and  are  solved  from  Equations  (3)  and  (6).  Substituting  Equation 
(7)  into  Equation  (3),  R"(0)  *  1, 


°  M2 


(8) 


Equation  (7)  can  be  written  as 

R(x)  -  <*Q  (x4  -  x3  (2M  +  l)  +  x2  (M2  +  2M)  -  xM2] 

+  [x5  "  x4  (2M  «■  1)  +  x3  (M2  +  2M)  -  x2M2] 


(9) 


Using  this  form  of  R(x)  and  applying  Equation  (6), 
1 


/  ROO  dx  ■  •„  [I  - 

1  (_6  5 


1  9 

<2M  +  1)  +  ~  (M2  +  2M) 
3 


1  9 

(2M  +  l)  +  -  (M2  +  2M) 
4 


2  *] 

r"2]-0 


Then,  applying  Equation  (8)  and  solving  for  aj ,  . 
a  -10M  +  10M  -  3 


1 


M2  (-5M2  4-  6M  -  2) 


(10) 


(ID 
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EVALUATION  OF  P(x) 


(12) 

(13) 

(14) 

(15) 

EVALUATION  OF  Q(x) 

Since  Q(0)  -  Q"(0)  -  Q(l)  -  0,  and  P(x)  is  a  polynomial  of  the  fifth  degree. 


Equation  (16)  becomes 

Q(x)  -  x2(x-l)  (Y0  +  Yi*  +  Y2x2)  (16) 

The  terms  Yg,  Yp  and  Y2  are  solved  from  Equations  (4),  (5),  and  (6).  From 
Equation  (4),  Q(M)  ■  1,  and 

YQM2  (M  -  1)  +  Y1M3(M  -  1)  +  Y2M4  (M  -  1)  -  1  (17) 

Also,  from  Equation  (5),  Q*(M)  ■  0,  and 

YqM(3M  -  2)  +  YjM2  (4M  -  3)  +  Y2M3  (5M  -  4)  -  0  (18) 

Equation  (16)  can  be  written  as 

Q(x)  -  Y0  (x3  ~  x2)  +  Yi  (x4  -  x3)  +  Y2  (x5  -  x4)  (19) 


Since  P(0)  -  P'(0)  -  P(M)  -  Prf(M)  -  P(l)  -  0,  and  since  P(x)  is  a 
polynomial  of  the  fifth  degree,  Equation  (12)  takes  the  form 

P(x)  -  x2  (x  -  M)2  (x  -  1)  B0 

The  term  gQ  is  solved  from  Equation  (6).  By  expanding  Equation  (12), 
P(x)  -  B  [x3  -  x4  ( 2M  +  1)  +  x3  (M2  +  2M)  -  x2M2] 


substituting  into  Equation  (6), 

1 

,2 

j  r\x;ax  —  v 

0 


f  P(x)dx  -  C' 

-  B0  [-  -  -  (2M  +  1)  +  i  (M2  +  2M)  -  -M2  1 
0  L6  5  4  3  J 


Solving  for  B0, 


B0  -  C 


7[ir-  j 


(2M  +  1)  +  ~  (M2  +  2M) 
4 


-M 
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Finally,  from  Equation  (6),  /  Q(x)dx  ■  0,  and  Equation  (19), 

0 

A’04*  ■  (r -  ?) + Y>  (r -  r) +  Tj ( * ' 0 ‘ 0 


From  Equations  (17),  (18),  and  (20),  the  system  of  three  equations  and  three 
unknowns  can  be  simplified  and  solved  simultaneously  for  yq,  Yj,  and  Y2*  The 


resulting  values  are 


1  -  2.5M2  12. 5M2  -  10M2  -  3M  +  2 

(1  -  1 . 5M)M  7.5M3  -  6M2  -  4M  +  3 


1  -  Y0  (1  -  2.5Mz)  (M  -1) 
M3  (M  -  1)  (1  -  1.5M) 


Y2  -  -  (2.5  Yq  +  1.5  Yt) 


The  set  of  equations  to  describe  an  axisymmetric  body  by  means  of  a  fifth- 
order  polynomial  can  be  summarized  as 

y2  -  2r0R(x)  +  CvP(x)  +  CZQ(x)  (24) 

where 

R(x)  ■  x(x  -  1)  (x  -  M)2  (aQ  4-  a^x) 

P(x)  -  x2(x  -  M)2  (x  -  1)  B0  (25) 

Q(x)  -  x2  (x  -  1)  (Y0  +  Yjx  +  Y2x2) 


°  M2 


.  -10M2  +  10M  -  3 

- - 

M2  (-5M2  +  6M  -  2) 
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1  -  Y„  (1  -  2.5M2)  (M  -  1)M2 
Yj  -  - 5 - 

M3(M  -  1)  (1  -  1.5M) 

\  -  -(2.5YC  +  1.5Yj) 

As  noted,  not  all  combinations  of  rQ,  C^,  C,  and  M  define  shapes  that  are 
either  realistic  or  smooth.  To  ensure  a  satisfactory  shape,  equations  have  been 
defined  that  must  be  satisfied.  These  restrictions,  which  are  the  same  as  those 
imposed  in  Reference  7,  include  the  existence  of  a  single  maxima  between  the 
nose  and  tail  of  the  body  and  at  a  prescribed  location,  M.  Also,  y  must  always 
be  greater  than  zero,  and  the  slope  of  the  body  must  continuously  be  decreasing 
as  x  goes  from  nose  to  tail. 
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